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ABSTRACT 

The    distribution   of   inventory  deaiand  over   procurement 
leadtime   is   modeled   using  mixed    probability   distributions 
that    explicitly   account   for    the    high    incidence   of   zero 
demands   observed   in   Inventory  Control    Point   Demand   History 
Files.      Analysis   was   limited   xo    the   right-hand   tail   area   of 
the    distribution   on   the   assumption   that    that   area    is    of   crit- 
ical   importance    in   determining   the   effectiveness    of   an    inven- 
tory  system.      Probabilistic    models   studied   were:       1) 
Bernoulli-exponential r   2)    B ernoulli-lognormal,    and   3) 
Bernoulli-logistic.      These    compound    distributions    were   com- 
pared  to   several   standard   distributions    including   the   Pois- 
son ,    negative   binomial  and    normal   distributions   using   a 
resampling    procedure    appropriate   in    ::ases   such  as    this    where 
a  paucity   of   data   exists.      Fits   obtained   from   the   75th   to 
95th    percentiles   indicated    the   mixed   distributions    may   be 
superior   as   a   class   to  the    standard   distributions    for   high- 
demand   items. 
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I.   INTRODUCTION 
Inventory  system  performance^  in  a  situation  with  random 
demand  depends  upon  reorder  points  computations.   Reorder 
point  computations  depend  upon  the  probabilistic  model  chosen 
to  represent  the  inventory  system  in  that  the  reorder  point 
is  composed  of  two  parts:   1)  The  expected  demand  during  a 
procurement  leadtime  and  2)  The  safety  level^  determined  from 
the  probabilistic  model  of  demand  over  leadtime.   Probabilis- 
tic models  are  generally  only  utilized  to  represent  demand 
over  procurement  leadtime  since  that  is  the  only  time  period 
when  stockouts*  potentially  occur.   The  chosen  distribution 
is  then  a  conditional  distribution  given  procurement  lead- 
time.   The  Navy  has  historically  used  three  distributions  of 
demand  for  models  utilized  at  the  Inventory  Control  Point 
(ICP)  echelon  of  the  Navy  Supply  System.   The  current  proba- 
bilistic models  and  the  average  annual  demand  that  is  used  to 
determine  which  model  to  utilize  are  displayed  in  Table  l.s 


1  See  Appendix  A  for  definition  of  "System  Performance". 

2See  Appendix  A  for  definition  of  "Reorder  Point". 

3S€e  Appendix  A  for  definition  of  "Safety  Level". 

♦See  Appendix  A  for  definition  of  "Stockouts". 

SThe  Aviation  Supply  Office  (ASO)  has  recently  changed  to 

the  use  of  the  normal  distribution  as  the  model  of  choice  for 


Table    1:      Current    Inventory   Control   Point 

Probabilistic   Deoaand  Over   Leadtime   Models 

Average   Annual 
Distribution  Demand_Ran5e 

Poisson  0-1         Low 

Negative  Binomial         1-20      Medium 

Normal  >    20      High 

If   the   actual   quarterly    demand   records    for  an   ICP   are 
examined,    one   striking  characteristic    that    is   present   in 
nearly  all    records    is   the  high   incidence   of   zero    observa- 
tions.     This   is   not   surprising   when   one    considers    the   ech- 
eloning  of    the   supply    system   and  the   fact   that   users   only 
place  a    demand  on   the    wholesale   system    as   their   reorder 
points  are    reached-       Since    most    activities   stock    at   least   one 
quarter's   expected    demand   for  an   item,    reorders   are   expected 
at   most   four   times    a   year.       By   the   natural    phase    differences 
of   each   activity's    reorder    actions,    the    ICP   may   experience 
zero    demands   during   a    giver\    quarter   for    even    high-demand 
items,    and   will   surely  often   experience    quarters    of    zero 
demands   for   medium-   and   low-demand   items.      A   previous   study 
of   Pacific   Fleet   Combat    Stores  Ships'    demand   data   [Ref.    1] 
found   that    the   frequency   of    zero    observations   in    those   demand 


all    demand    categories.      Ships  Parts   Control   Center    (SPCC) 
still    uses    all   the    distributions    in    Table    1. 
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records   appeared  to   have   a    low  correlation   with   the   overall 
demand   level    (i.e.    high,    medium    or   low).* 

A   desirable   model    for   an  inventory   system   should    compute 
reorder   levels  that   accurately  correspond   to   a  specified 
stockcut   risk   based   on  the    level   of   actual   or   anticipated 
demand   for    each    item.      If  such    a  model    were   available,    the 
inventory  system   could  avoid   the   added    expense  of    either   hav- 
ing   too   much   or    too   little    stock.      One    problem   with   the   dis- 
tributions   in  Table    1    is  that  none   can   account   for   the 
probability  mass   at   zero;    therefore,    none   can  accurately   com- 
pute  reorder   levels. 

This  study   explores   the    use    of   compound   probability   mod- 
els   and   their   potential   to    more    accurately   achieve   the   goals 
presented  above    by   explicitly  accounting   for. the    probability 
mass    at    zero.      The   results    obtained   should    be   considered   as   a 
first   step    in   the   exploration   of   a   class    of   models   which   have 
heretofore    received   little    attention   in   the   context    of    model- 
ing   inventory   demand. 


*See   Section   II. B.I    for    a   more   detailed   discussion   of 
this    study. 
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II .       BACK GROJND 

A.       CURRENT    MODEL   SELECTION    TECHNIQUES 

In  the    past,    the    usual    approach   ot   choosing   the    probabi- 
listic  model    for   lead   time    demand   was    based   on   matching  the 
empirical  cumulative   distribution   function   determined   from 
the    actual   demand   to    various   theoretical   cumulative   distribu- 
tion   functions.      The    shapes    of   the   two   curves   are   compared 
using  the   Kolmogorov-Smirno v   Statistical   Goodness    of   Fit 
test.''^   The   Chi-Squared  test    has   historically  not    been   used 
because    it    is  not   very  powerful    when   the   number   of   observa- 
tions is  small    {supply  demand   data   is   retained   for   at   most 
twelve    quarters    and    usually    only    eight    quarters    at    the    ICP 
level).      The   Kolmogorov-Smi rnov    test    provides   only   a    relative 
measure   of    goodness   of  fit,    but    it   is    usable   with    very   small 
sample   sizes.      The   goal   then   is   to   select   the   model   with   the 
best   agreement   between  the    empirical   and   theoretical   cumula- 
tive   distribution   function    curves- 

The    Navy   Fleet    Material    Support    Office    (FMSO)     has    uti- 
lized  the   above    Kolmogorov- Smirnov   test    to    evaluate   several 
theoretical  distributions.       Most   recently   [Ref.    2]   FMSO 


^A   good   description   of    the   Kolmogorov-Smirnov   one   and 
two-sample   tests   may    be    found   in   Siegel,    S.,    Non  _  Parametric 
Statistics,    p. 47-52    and    127-136,    McGraw    Hill,    1956. 
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evaluated  seven  distr ibutionss  at  the  90"^  confidence  level 
using  demand  data  from  NAS  Brunswick.   The  use  of  demand  data 
from  a  stock  point  vice  an  ICP  assumes  the  underlying  demand 
behavior  at  the  stock  point  and  ICP  are  similar.   Though  this 
was  never  validated  directly,  the  NAS  Brunswick  demand  data 
were  deemed  to  have  "unique  qualities  better  suited  for 
demand  analysis"  as  follows: 

1.  Twelve  quarters  of  da mand  and  demand  frequency  data 
were  available  per  item,  SO'^  more  than  available  in  ICP 
files, 

2.  New  item  identifiers  that  permit  the  selection  of  only 
steady  state  demand  items  were  available, 

3.  Reliability  of  the  data  was  established  from  prior  FM30 
studies. 

The  FMSO  analysis  showed  that  Navy  demand  patterns  are 
very  poorly  modeled  by  any  of  the  seven  distributions  stud- 
ied.  No  new  models  were  proposed  and  no  attributes  of  the 
demand  patterns  were  indicated  as  the  most  significant  factor 
in  the  failure  of  the  standard  distributions  to  model  Navy 
demand  patterns. 

There  are  several  problems  wi-h  the  way  the  Kolmogorov- 
Smirncv  test  was  used  in  the  FMSO  study: 


sDistr ibutions  tested  were:  Normal,  Negative-Binomial, 
Poisson,  Logistic,  LaPlace,  Gamma  and  Uniform. 
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1.  The  procedure  is  not  strictly  applicable  in  the  case  of 
discrete  distributions  such  as  the  Poisson  and  negative 
binomial. 

2.  The   test   assumes  that    the    parameters   of   the    distribu- 
tion   being  tested  are    completely   specified    in   advance, 
i.e.    not   estimated   from   the   data.       (Lilliefors    [Ref.    3] 
developed   alternate   tables    for   use   in   testing   the   expo- 
nential distribution    when    parameters   are   estimated    from 
the  data.      This    procedure    was   used   in   Reference    1.) 

3.  The  Kolmogorov-Smirno V   test   evaluates   goodness-of-f it 
over   the   entire   range    of  the   distribution.      For    inven- 
tory   problems,    the   region    of   greatest   interest    is   the 
right-hand   tail   of   the   distribution. 

B.       PREVIOUS    RESEARCH 

Research   into   distributions   which    have    properties   that 
more    closely   model    actual   demand   patterns   has   been   very   lim- 
ited  and,    when  conducted,    often   did   not    explore    anything   but 
the    usual   standard   distributions.      Two    projects   though   stand 
out    because  of  their   innovative    use    of    non-standard   distribu- 
tions and  the   good    fit  achieved. 

1 .      Bernoulli-exponenti al  D  istribution 

In    Reference    1,    Pacific    Fleet    Combat   Scores    Ship 
demand   was    modeled   using   the   mixed   Bernoulli-exponential 
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distribution.      This  distribution    was   chosen    because    it    was 
observed  that   nearly    all    line  items   had   an    unusually    high 
occurence  of   demands   of   zero.      Analysis   of   the  data   indicated 
that    with   high    probability,    the    number    of   non-zero    observa- 
tions   for   a   given    item   was    unrelated    to   the   average   value    of 
those   non-zero   observations.      The   conclusion   was    that    the 
demand   process   could    be  thought    of   as    two   independent    subproc 
esses   with    one   process  determining   whether   a   demand   would 
occur   or  not   and  the    other    determining    the   quantity   of   the 
demand,    given   the   demand   did   occur.      The   former   process    was 
modeled    as    a   Bernoulli  process   with    parameter   p,    p    being  the 
probability  that   a   demand  did  occur.      The   latter    process   was 
modeled   as    a   continuous  exponential   process   based    on    explora- 
tory  data  analysis.      The   resulting   probability  distribution 
function   for   the   number   of    units   demanded   in   a   leadtime    for 
the    Bernoulli -exponential  distribution    is   given    by: 


1    -    p 


,    X    =    0 


(1) 


h(x)    =   < 


pXex  p  (-Ax)  ,    X    >    0, 


where : 


VX  = 


P  = 


the   expected   value   of    demand,    given   that   the 
demand    is   greater   than   zero, 
probability  that    demand   does    occur. 
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Th9  complementary  cumulative  distribution  function  is  given 
by: 

H  (X)  =  p  exp(-Ax)  ,   X  >  0,  (2) 

where  H  (x)  is  the  probability  that  demand  will  not  exceed  x; 
this  is  equivalent  to  the  complementary  cumulative  distribu- 
tion   function   used   by    Hadlay  and   Whitin   [  Ref .    4]. 

The   hypothesis  that    the    demand    of   Combar    Stores   Ships 
could   be   modeled  according    to   the   above    probability   distribu- 
tion   was  tested    using    the   K olmogorov-Smirnov   test    described 
earlier   plus  tests   on    the  theoretical    risk   and  the   observed 
risk    for   values   of   risk   in    the   upper   right   hand   tail    area   of 
the   distribution.      The  results   of   these    tests   on    five    inde- 
pendent  samples    from   the   demand    data   available  provided 
strong   evidence   that    the   Bernoulli-exponential   distribution 
describes   the   demand    for   any   given   stock   item   very   well. 
2.      loqar  ithmic-Poisson-Garama   Distribution 

A  model  for  the  distribution  of  demand  during  pro- 
curement leadtime  [Ref.  5]  was  derived  under  the  following 
assumptions: 

a.  Requisitions    occur   according   to    a   stationary   Poisson 
Process, 

b.  Requisition   sizes   follow  a   logarithmic   distribution. 
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c.      Leadtime   is   a   random    variable   with   the   gamma 
distribution. 
The   rssalting   probability  function   for   the   number    of   units 
demanded   in   a   leadtime   is: 

X  ^  k 

k  r— »  (3) 


""  -(:A--^^--^.  n"-^-i)^>=.^D^- 


D 
k=l  3=1 

for    X    =    1,    2t    3,    ... 

and 

where: 

a^6  =  parameters  of  the  Gaami  distribution  defining 
leadtime, 
0  =  parameter  of  the  logarithmic  distribution 
solved  for  by  interval  bisection  from 
E  (X)  =  -0  /  (1  -  0)la(1  -  0)  , 
X  =  requisition  arrival  rate, 
c  =  -A/  ln(1  -0) 
^x  k ~  Stirling  numbers  of  the  first  kind  computed 
from  the  recursion 


^x,k~  ^x 


-l,k-l  -^  (^  -  '')Sx-i,k 


for  k  =  1,  2,  ...,  X   and  x  =  1,  2,  ... 

with  Sj^  Q  =3  for  all  x 
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This  distribution  is  called  the  Logarithmic-Poisson-Ganinia 
(LPG) .  It  is  defined  from  four  parameters  a,  6,  Q,  and  X 
(0and    X  determine   c)  . 

For  more   efficient    computations,    a   recursion    equation 

for   integer    a  was   developed. 

X 

h(x)    =   Y^'^x^k  ^"^ 

k=i 
where: 

^x,k    =    'r   ^^^2^°'   *    ^    '    ^)^x-l,k-i-^  (^    -    ^)^x-l,k3 
C2    =   c   /    (X  +  6) 

and    0,    c,    a,    6   and    X  defined  as    in     (3). 

Computing  reorder  points    from    fractiles   of   the   LPG 
distribution,    even    using   the   simplified    version   of    (4) ,    is 
generally   too   time   consuming   for    many    real    applications.      As 
a  result,    an   approximation    was  developed   and  tested   that    used 
a  scaled  version   of   the   Poisscn    distribution   to   approximate 
the    negative   binomial    distribution.      The    resulting   approxima- 
tion   is : 

(a   *    X    -    1)  !     /      s      \7     M      X""  (5) 

p{Z(t)    =   kx}     =      —7-. 

x!  (a  -     1)!         \li   +    S/\U   +   6/ 

for    X    =    0,    1,    2,    ... 

where: 

Z(t)    =    number   of    units    demanded   in    time   t 
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M  =  C0 

a,  6,  c  and  0  defined  as  in  (3). 

The  authors  conducted  subjective  tests  of  goodness  of 
fit  to  their  LPG  model  and  its  derivatives  using  Air  Force 
consumable  items  and  concluded,  using  as  an  example  one  item, 
that  there  was  a  "very  close  agreement  between  the  observed 
and  the  predicted  cumulative  distribution  functions  for  this 
item.  " 
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III.       PROPOSED_MODELS 

A.       COMPOUND    BESNOOLLI    MODELS    IN    GENERAL 

The   results    of   the  previous   research   indicate   that   a 
potential   exists   for   improving   inventory   performance   by  the 
use    of  slightly    more   complex   models.      The   focus   of   this   study 
is   on   the  compound   Bernoulli   models,    leaving   further    analysis 
of   the  LPG    model  and    its   derivatives   to    others.      The   use   of 
compound   Bernoulli    models  has  several    advantages    over   stan- 
dard  simple   distributions: 

1.  The   model,    by   necessity,    has   an    added   parameter   whose 
sole    function   is  to   estimate  the    probability   of   a    non- 
zero   demand, 

2.  For   simple   models   defined    only    over   the   positive    real 
numbers,    the   addition    of   the   Bernoulli   parameter   simply 
results   in   the   scaling  of    the   cumulative  distribution, 
thus    adding   little   additional  complexity   to    the    model. 

3.  By  explicitly   accounting   for   zaro    demands,    compound 
models   account    for   more  of   the   observed   variance   than 
do    simple    models. 

Compound   models   are   not    without   disadvantages: 
1,      Each    parameter   added    to  the    model    must    be   estimated   for 
each    individual   stock    item.      For    a   large  inventory 
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system   such   as   the   Navy's,    this   equates   to    considerable 
additional  computational   requirements. 
2.      Additional   parameters    must   be    stored    for   future 
retrieval    or    else   computed   from    stored   data    when 
required,    thus   either    extending   computational   time   or 
else    requiring    additional   online   storage   space. 
The   Bernoulli   compound    models  all    have   similarly   shaped 
cumulative    distribution   functions.      Figure    1    is   a    sample   dis- 
tribution  function    for  the    compound   Bernoulli-lognormal   dis- 
tribution  with   a   Bernoulli    factor   of    p    equal   to    0.6,    a    mean 
of    1.0   and   a   standard    deviation   of   0.5.       Note  that   the   dis- 
tribution  has   a    mass   of   zero   for   negative   observations,    a 
mass   of    1-p   at    zero   and   the    usual   cumulative   distribution 
function   shape   for    positive    observations   but   of   mass    p    vice 
mass    one. 

B.       MODEL    DESCRIPTIONS 

Three  compound   Bernoulli   models   were   formulated    for   eval- 
uation.     Each   model   was   derived    from    the   base   distribution 
after   compounding    with  the    Bernoulli    process.      kn    additional 
requirement    was    made    that   the  inverse    of    the   compound    model 
must    be   computationally   relatively   simple. 

1  •      Bernoulli-exponenti al  D  istribution 

This   distribution   was  derived    in    Reference   1    and  the 
derivation    will    not   be  repeated    here.       Equations     (1)     and    (2) 
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Figure    1:      Compound   Cumulative   Distribution   Function, 
Bernoulli-lognormal   Distribution 

are   the   probability   density    and   complementary  cumulative 
distribution   functions  respectively.      The   Bernoulli- 
exponential   inverse   complementary   cumulative   distribution 
function   is   computed   as: 


X    =  < 


r      H(X)      >     p 


M[ln(p)    -   In  (H(x))  ],    H(x)     <   p 


(6) 


The    parameters   of  this  distribution   are    estimated   from   the 
sample   data   as   follows: 
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p  =   number    of   non-zero    demand   observations    in   sample 

divided   by   sample    size, 
y  =   X   /  p,    where   x   is    the  entire    sample   mean. 


or 


y   =   X,    where  i   is  the    sample   mean    of   only   the  non-zero 
demand   observations. 
2 .      Bernoulli -log normal    Distribut ion 

The  lognormal   distribution   was   chosen    as    a   candidate 
for   testing   because   of  several   desirable   properties: 

a.  It   is   defined   for   all    real   numbers    greater    than   zero, 

b.  Its   inverse    is   readily   computed    using   inverse   normal 
approximations    that   are  well   documented   for    accuracy, 

c.  Alone,    the   lognormal    distribution    does  not    allow   the 
case    of   zero    demand,    but   when   combined   with    the   Ber- 
noulli  distribution,    the   compound    Bernoulli-lognormal 
distribution    defines    all   the   expected   demand   values 
greater   than    or   equal    to   zero. 

Utilizing   the    Bernoulli    parameter   p   as   a    scaling   fac- 
tor   for   the   lognormal   density  function,    the   compound 
Bernoulli-lognormal   density    function    may    be    expressed   as: 


1    -   P 


,    X    =   0 


(7) 


h(x)    =    { 


/ln(x)     -  y\ 
P<t> U    X  >   0 


where : 
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p  =   probability    of   a  non-zero   demand, 

y  =   mean   of   the    natural   log   of    all    demands   greater   than 

zero, 
a   =   standard   deviation    of   the   natural    log   of    all   demands 

greater   than   zero, 
<t)   -   Standard  normal   density   function. 
Integrating    (7)    results   in    the  complamentary   cumulative   dis- 
tribution  function   given   by: 


,    X   =    0 


(8) 


H(x)    = 


In  (X)    -  yvi 


1    -    $ 


X    >   0 


where: 

$    =   Standard   normal  cumulative   distribution    function. 
The    Bernoulli-lognormal   inverse    cumulative    distribution   func- 
tion   is   easily   derived  from     (8)     as: 


X    =< 


r     H(X)      >      p 


exp[y    ^     a$-i  {1    -    H  (x) /p)  ],    H  (x)    <    p 


(9) 


The    parameters   are    estimated   from   the    sample   data    in    the 
usual   manner:  . 

p  =   number   of   non-zero   demand   observations   in   sample 
divided   by   sample    size. 
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y   =   X,    where  x    is  the    sample   mean    of   the   natural   log   of 

all   non-zero   demand  observations, 
o  =  s,    where  s   is  the    sample   standard   deviarion   of  the 
natural   log   of   all    non-zero    demand  observations.      An 
alternate    method  for  estimating   the   mean    and   stan- 
dard  deviation   of    the   lognormal   distribution   is   to 
use   the   method   of    moments.      The   procedure   in    effect, 
transforms   the    mean  and    standard   deviation   obtained 
from   the    unt ransf o rmed    sample    data   vice    transforming 
the   data    first   and    then   computing   the  sample   mean 
and   standard    deviation.      This    method   involves    fewer 
logarithms   and   will   run    faster   on   a   computer.      The 
procedure   is    described    in   Appendix   D. 
3.      Bernoulli-lcqistic    Distribution 

The   logistic   distribution   is   a    pseudo-normal    distri- 
bution  that    is   similar  in   shape   to  the   normal,    but   is    easier 
to   handle  mathematically.      Its   density    function    is   defined   as 
[  Re  f .    6  ] : 


f(x)    = sech^i 

4\/3a 


7r(X    '   uY 

^     2yTa   , 


,    -«      <    X   < 


(10) 


while    its   cumulative   distribution   function    is   defined   as: 


1  r                     /7r(x   -   y  ) 
F(x)    =   -    1    -I-    tanh 

2  1  \      2\J3  0 


(11) 


,    -«     <    X    < 
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where: 

\ji  -   mean  of   the   en-cire    sample, 
a  -  standard   deviation    of   the   entire   sample. 
With   the   logistic   distribution,    it   is    no   longer   a   simple   mat- 
ter  of  scaling  the   distribation    by   the    Bernoulli    p   parameter, 
since   this    distribution   is    defined  over    the   entire   real    line 
and   not    just   on   the    positive   half.      Consequently,    the   distri- 
bution  requires    a   new    constant   of   integration   to    replace   the 
1/2    used  in    (11).      Integrating    (10)     from    0   to    »    provides   the 
new   consxant   as: 


/7r(x    -  u)\         "                      /  ^y  \ 
sech2( ;;:: )dx   =    1    •»■   tanhf 


(12) 


which,    when   compounded   with    the    Bernoulli    distribution,    gives 
the    Bernoulli-logistic  distribution   density   function    as: 


/  1   -    p 


h(x)    = 


/   Try    > 

1    >    tanhl 

\2\!3o^ 


sech2| 


•7r(X    -   M  )■ 

-      2^a   ^ 


,     X    =    0  (13) 


,    X   >    0 


,    otherwise. 


where : 


u  =   mean   of   all   non-zero  demand   observations. 


26 


a  =   standard   deviatiori    of   all   nDn-zero   demand 
observations. 
After   integrating    (13),    the    complementary    cumulative    distri- 


bution function   may   be, written  as: 

r 

1  r       /  Ty  \  /^(x  -  y  )• 

H(x)    =   p{^    -     )tanh— --  U    tanh z:^— 

1    ♦   tanhf \ 


\2/ra/P 


(I'*) 


for   X   >    0 


The    inverse  cumulative  distribution   is: 


-{ 


/so 
u  ■•• In^ 


2p 


V. 


H(x)     1    *    tanhf—--    , 


,    H(x)     >    p       (15) 


-     1    >    r     H(X)       <     p 


The   parameters   are   estimated    from    ths    sample   data    in    the 
usual   manner: 

p  =   number    of   non-zero    demand   observations    in   sample 

divided   by   sampla    size, 
y  =   X,    where   x    is   the    sample   Dean    of    all   non-zero   demand 

observations , 
a  =   s,    where   s    is   the    sample   standard    deviation   of   all 
non-zero   demand   ob servations. 
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IV.       EVALUATION_PHOCEDURE 

A.       THE    DATA 

In   order  to    test    the   validity   of   the   proposed    models, 
samples   of    actual   demand   ware   obtained    from   the   Operations 
Analysis   Department   at  FMSO.      The   data,    accumulated    from   the 
demand  history    files    of   the    Aviation    Supply   Office,    was   orig- 
inally  used   as   input   data   for  the  5A    (Aviation  Afloat   and 
Ashore  Allowance   Analyzer)     [Ref.    7].      The   data   consists   of 
1587    consumable    1R    cog   items   and    2892    non-program-related' 
repairable    2R  cog    items.      The  information    for    each   item   is 
contained   in   one   master    record  and   several   subrecords.      The 
master   record   contains   identifying   information   on    each   item 
such    as   the   national   stock    number,    replacement   price,    etc. 
Each    subrecord   contains   up    to   forty-six   demand   records,    each 
including  the   demand    quantity  and   the    day   of  the    demand.      A 
complete  record   layout   is  contained   in    Appendix   E.      To   ease 
processing,    the   Julian   Dates   of   the   original   demands    were 
replaced  by   FMSO    with   a    sequential   date   ranging   from    1    to 
1500    representing   approximately    four   years    of   available 
history. 


'See   Appendix   A   for   definition   of    "program  related" 
items. 
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For   this   application,    additional   screening   was  conducted 
to   group  the  demands   into  thirty-day    buckets,    thus    providing 
a  demand  time   series    of    forty-eight   observations.      Further 
editing    was   required    prior    to  analysis   to   remove   the    negative 
demands  noted   for    many  of   the  items.      It   was   assumed   that    the 
negative  observations    were    the   result    of   cancellations    of 
previous  demands   which  were    never   filled.      Therefore,    when- 
ever   a  negative   demand  was    encountered,    the   demand   series    was 
searched  backwards    for  the    first   positive   demand    equal   to   or 
larger  than   the   negative   demand.      If    such    a   demand   was    found, 
the    positive   demand   quantity   was   reduced   by   the   absolute 
value   of  the   negative    quantity   and   the    negative   quantity    was 
set    equal  to   zero.      If  no   offsetting    positive   demand    was 
found,    the   negative   demand    quantity   ^as    still  set    equal   to 
zero.      The    resulting    edited    demand   time    series   was   used   as 
the    base   for   all  further   analysis.      Two    additional   demand 
series   were  created   from   the   edited   series: 

1.  k    demand  series    of   only   the   non-zero    demands, 

2.  A    demand   series    of   the   log   of   the   non-zero    demands. 

The  sample  mean  and  sampla  standard  deviation  for  each  of  the 
three  demand  series  were  oomputed^o  ^^^  the  group,  represent- 
ing   one   line   item,    was  accepted    for   futher   processing   if: 


1 ostationarity    was   assumed   for   the    demand  time   series 
under   study.      If   a    trend   ware   in    fact    present,    the   resampling 
technique  effectively    would    eliminate    it    by   shuffling   the 
observed  demands   into    random   order. 
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1,  There  were   at   least  two  non-zero    demands,    and 

2.  The   standard   deviation  of    all   three   series    was   non- 
zero. 

Once  accepted,    the  sample   mean   and/or   sample    standard 
deviation   of   the  appropriate   series   was    used   to   estimate   the 
value   of   X,    such   that    for   a    given   probability,    demand    will 
not    exceed   that    value   of   x.       The   values   of    x   were   computed 
for    each   model   tested    at   several   probabilities   from   the 
inverse   cumulative    probability  functions.      Table    2   displays 
the    sample    mean   and/or  sample  standard    deviation    required   for 
each    model    investigated. 

B.       THE    RESAMPLING    PROCEDURE 

The   resampling   procedure   used   was    chosen   over   the   more 
traditional   methods   discussed   in   Chapter   II   because    it    pro- 
vides  a    method   of   comparing    a  theotetical   distribution   to   the 
sample   data   at   specified    probabilities    and   is  applicable   in 
cases   such   as   this    where  there  is   a   paucity   of   data   available 
for    analysis.      This   procedure   is   similar   to   the   "bootsxrap" 
procedure   used   by   Efron   [Ref.    8].      The    idea   behind   zhis   pro- 
cedure  is  tc   randomly   sample    with   replacement   the    series    of 
available   data   to   create   additional   pseudo-samples   that   pos- 
sess  the  same  statistical  properties   as    the   original    sample. 
The    desired   statistical    property   can   then    be   estimated   from 
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Table  2:   Required  Data  for  Inverse  Computations 


STANDARD 

MODEL 

MEAN 

DEVIATION 

1- 

Exponential 

A 

N/R 

2. 

Normal 

A 

A 

3, 

Poisson 

A 

N/R 

4. 

Negative   Binomial 

A 

A 

5. 

Lognormal 

C 

C 

6. 

Logistic 

A 

A 

7. 

La  Place 

A 

A 

8. 

Bernoulli-exponential 

B 

N/R 

9. 

Bernoulli- lognorma 1 

C 

C 

10, 

Bernoulli- logistic 

B 

B 

Where: 


A    is   from    edited   demand   series 

B    is    from   non-zero    demand   series 

C    is   from    log  of  ion-zero   demand   series 

N/R   indicates   parameter   Not   Required   for    this    model 


each    pseudo   sample.       Wirh  a    sufficient    number  of    repetitions, 
the    distribution   of   the    estimated  property   is   known    to    be 
normal   from   the    central   limit    theorem    with   mean    equal  to   the 
population    mean    and   standard   deviation    equal   to   the   popula- 
tion   standard   deviation   divided    by   the    number  of    repetitions. 
Thus   by   the   use   of  the    resampling    procedure,    it.   is    possi- 
ble   to   study   theoretical  distributions    at   points    in   their 
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right  hand  tail  area,  which  is  the  region  of  most  interest  in 
computing  safety  levels  in  inventory  models.   The  resampling 
procedure,  as  applied  in  tais  study,  is  a  relatively  straight- 
forward application  of  ths  following  technique: 

1.  The  edited  series  for  each  sample  of  n  demand  observa- 
tions is  treated  as  the  population  from  which  the  ran- 
dom samples  are  drawn. 

2.  A  random  sample  of  size  n  is  drawn  from  this  popula- 
tion, creating  a  pseudo  sample.   The  pseudo  sample  is 
not  a  permutation  of  the  origiial  population  since  the 
sample  values  are  selected  with  replacement  from  the 
original  population.   This  is  easily  accomplished  in  a 
computer  by  generating  n  uniform  random  integers  in  the 
range  from  1  to  n,  and  using  these  numbers  as  sub- 
scripts to  select  the  sample  from  the  original 
population. 

3.  The  pseudo  sample  percentile  is  computed  for  each  rep- 
etition from  (16)  and  is  compared  to  the  theoretical 
percentile. 

p  =  (number  Df  demand  observations  <  x)  /  n    (16) 

The  various  values  of  x  are  computed  from  the  inverse 
cumulative  probability  functions  for  each  model 
evaluated. 
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Dnder  the  null  hypothesis  that  the  sample  demands  come 
from  a  particular  probability  distribution,  the  expected 
value  of  p  should  equal  p,  the  theoretical  percentile: 

EC  p]  =  0  under  H^  (17) 

and   the   distribution   of    p  should   be  Normal  (p,a2/n) .      The   per- 
centile   p  is   estimated  as  above    and   a^    is   estimated    from 
standard  binomial   results   as: 

VarCS]  =   p(1    -  p)    /  n.  (18) 


C.   MEASURES  OF  EFFECTIVENESS 

When  the  above  procedure  is  executed  at  several  theoreti- 
cal percentiles,  the  fit  of  the  sample  data  may  be  compared 
at  several  locations  in  the  theoretical  distribution  and  the 
overall  fit  evaluated.   Several  common  measures  of  effective- 
ness are  available  to  provide  a  criteria  on  which  to  judge 
the  success  or  failure  of  a  model.   Some  of  those  available 
are : 

1.   Maximum  absolute  deviation, 

2-   Mean  squared  error, 

3.   Algebraic  sum  of  errors. 
Method  three  was  eliminated  since  errors  may  offset  one 
another  resulting  in  a  seemingly  good  fit,  but  in  reality  a 
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very    poor   fit.      Method  one    is  potentially   a   good    way    to   eval- 
uate  error,    but    is   generally   not   as   sensitive   as    method   two. 
Therefore,    the   mean   squared    error  of    the   pseudo   sample   per- 
centile   estimates   was   chosen  as   the   measure  of  effectiveness 
of 'Choice   and  wcls  accumulated  as: 

^^  (19) 


i=l 


D.   COMPUTER  VALIDATION 

The  computer  programs  for  this  analysis  were  written  in 
FORTRAN  IV  for  execution  on  the  IBM  3033  attached  processor 
(System  370)  computer  installed  at  the  Naval  Postgraduate 
School.   The  programs  are  batch  oriented  since  the  volume  of 
data  available  prevented  storage  directly  on  the  user's  pri- 
vate disk  space.   The  data  tapes  received  from  FMSO  were 
loaded  onto  the  system  mass  storage  device  and  when  required, 
all  or  portions  of  the  data  were  transferred  to  the  Virtual 
Machine  (VM)  facility  for  processing.   Two  programs  were 
written,  the  first  to  copy  the  data  to  make  it  accessible  to 
the  VM  system  and  the  second  to  conduct  the  actual  analysis. 
A  complete  listing  of  the  FORTRAN  source  code  is  available  in 
Appendix  F.   Each  program  or  subprogram  contains  a  documenta- 
tion block  at  the  top  that  provides  a  complete  description  of 
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the    program's   purpose,   a   definition   of    variables    used    and   a 
list   of   user-written   subroutines   and    functions  required. 

The   random   number    generator    ussd   is   a    part   of   the 
LLHANDOd   II    series    [Ref.    9]    developed    at    the   Naval   Postgradu- 
ate  School.      The   programs  and   subprograms    were   thoroughly 
tested  both   independently   and  as    a    unit    using   known    data. 
The    results   were   verified  by   hand   held   calculator   and    in   all 
cases   agreed   with   the    program  output. 
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V.        RESULTS 

The  results  of  the  distributional  analysis  conducted  on 
the  ten  distribution  models  listed  in  Table  2  are  tabulated 
in  Appendix  B  for  1R  consumable  items.  Appendix  B  is  divided 
into  the  three  demand  classes  listed  in  Table  1,  and  then 
each  class  is  divided  by  distribution  and  percentile  within 
each  distribution.  The  results  for  the  2R  repairable  items 
are    presented   in   a    similar    format   in    Appendix   C. 

The   consistency  of  the    resampling    method   with    varying 
data    was   evaluated   by   testing  independent   sections   of   the 
consumable    and   repairable   data   sets   and    comparing    the 
results.      In   both   cases,    independent   subsamples  were   created 
by   sampling   every   third   item    but   changing   the    initial   item 
sampled   from   the   first  to   the   third  item.      The   analysis,    not 
presented   here,    showed  that    the    resampling   method    produced 
consistent    results    ever   all    of   the  subsamples   tested.      The 
total   mean    square   error^i    of   the   distributions   in    each   sub- 
sample   varied   by   no    more   than   20%   from   the   total    mean   square 
error   values   obtained    from    the   entire    sample   population   as 
listed  in  Appendix    3    and   Appendix  C.      The   relative   standing 


1  iTotal    mean    square   error   is   the    sum   of   the    mean    square 
errors  computed    at    the  75th,    80th,    85th,    90th   and    95th 
percentiles. 
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of   each    distribution    was   not   altered    in   any   of  the 
subsamples. 

The   effect   of   the    number   of    pseudo    sample   repetitions    was 
also    studied.      The   results    proved   similar   while   varying  the 
number  of   repetitions   from    twenty  to    fifty.      The    standard 
deviation   of   the   percentile    estimates    decreased   as   the   number 
of   repetitions   increased.      The   total    mean   square    error    for 
each   trial   varied   randomly    but   remained    within   10^   of   the 
overall    values   for   each   distribution    with   more  repetitions 
producing  consistently  smaller   variations.      For    the   data 
analysis  production  runs,    forty    was   chosen   as   the    number   of 
resampling   repetitions  as  a    trade   off    between   computer   run 
time    and  minimizing   variance. 

A.       LOW-DEMAND    ITEMS 

A   summary   of   the   total    mean    squared    error   figures   for   the 
models   providing  the    best   fit   or   the   smallest  total   mean 
squared    error   is   provided   in   Table   3    for   the   low-demand 
items.      The   best    models    for    the    low   demand   items    all   had   very 
small   total   mean   squared   errors.      This   can   be   attributed   to 
the    very  low   probability  of    a   demand   ever    being   greater   than 
one    and   to    the   integer  sampling    plan    for   accumulating   the 
Bernoulli  trials.      For  the    integer   results    of   the    Poisson   and 
negative   binomial   distributions,    the   estimated   percentile 
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TABLE   3:        LOW-DEMAND    ITEMS, 

Total   Mean    Squared   Error    (lO-*) 

Negative 
Pol  sson        binomial        Comgoi^Ii^ 


Consumable    1R 

5.64 

7.86 

10.77 

Repairable   2R 

5.  86 

7.57 

12.91 

will   closely    match    the  theoretical   percentile   unless   the   num- 
ber  of   Bernoulli  successes    is   less   than   the   expected   number 
of   successes   for   a    given   theoretical   percentile.      This   pro- 
duces  a   conservative   result    that    will    tend   to   overstate  the 
quantity  of   stock   required    by  one   unit   to   provide    a  stated 
level   of   protection,    but   in    the   case   of   the   small    quantities 
involved,    is   not   an    undesirable   praparty.      The  compound   dis- 
tributions,   as   a   group,    gave  the   same   total    mean    squared 
error   for   much   the   same   reasons.      The    probability    of    a   non- 
zero   demand   is   very   small   for   the   low-demand   items,    generally 
less    than   0.05,    thus    these    models   tend    to    predict    a   quantity 
of   zero   for  all    percentiles.      The   resampling    procedure   and 
the   sampling   plan   again   provide    a   conservative  result. 

B.       MEDIUM-DEMAND    ITEMS 

A   summary  of   the   total    mean    squared   error   figures   for   the 
models   providing  the    best   fit  or   the    smallest   total   mean 
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squared   error   is   provided   in   Table  4    for   the    medium- demand 
items.      Here  again   the   predicted   demand    for   the   various   per- 
centiles  is   small.      The   resampling   technique   still   gives   the 


TABLE    4:       MEDIUM-DEMAND    ITEMS, 

Total    Mean    Squared   Error    (10-*) 

Negative  Bernoulli- 

binomial        Poisson        loqnormal 

Consumable    1R  28.55  48.28  80.60 

Repairable   2H  29.14  44.84  112.44 


edge    to   the   integer  distributions.      The    compound    distribu- 
tions  have    begun   to   show   some   differsnces   among   themselves, 
but    as   a  group   are   clearly    better   than    any    of   the   continuous 
simple   distributions.      The    magnitude   of   the   total    mean 
squared   error    is   significantly   larger    than    for   t.he    low-demand 
class    primarily   due   to  the    larger  allowable   total   demand 
range   defining   the   medium-demand   class. 

C-       HIGH-DEMAND    ITEMS 

The    summary   total    mean    squared   error    data    for    the    high- 
demand   items   are   listed   in    Table   5.       with    the   unlimited   range 
of   mean    demands   in   this   class,    the   higher   total   mean    squared 
errors  are    expected.      The   compound  models   have  come    into 


TABLE   5:       HIGH-DEMAND    ITEMS, 

Total   Mean    Squared   Error    (10-*) 

Bernoalli-        Bernoulli- 

lognor  mal  exponential        Exponential 

consumable    1R  168.25  188.90  298.21 

Repairable   2S  183.77  221.80  289.01 

their   own   in   this   class  with  the   Bernoulli-lognormal    yielding 
a  total    mean   squared   error    less   than    one   half  that   of   the 
normal   distribution.      The  second   best    distribution,    the 
Bernoulli-exponential,    gives  a   total   mean   squared   error   only 
slightly  higher   than   the   Bernoulli-lognormal. 
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VI.       CONCLUSIONS  . ANp„RECg MM ENDATIONS 

A.       SUMMARY 

This  distributional  analysis,    conducted   using   a   com- 
pletely   different   experimental   procedure   than  prior   analyses, 
provides  satisfyingly    similar  results    for   the  low   and   medium 
demand   categories    in   showing   that   the    Poisson   distribution 
gives  the   best   fit    for  the    low-demand    category   of    items   and 
that    the  negative  binomial    distribution    gives   the    best    fit 
for    the    medium-demand    category   of   items.      The   class   of   com- 
pound  distributions,    as   a   whole,    gave    good    fits    in   the   low- 
and    medium-demand   categories,    though   the   total   mean    squared 
error  was   two  to   four    times    as  large    as    that   of   the    Poisson 
or   negative   binomial   distribution. 

The    analysis   indicated    that    there    are   several    distribu- 
tions   which   provide   better    results   than    the   normal   distribu- 
tion   for   high-demand   items.       The    Bernoulli-lognormal 
distribution   consistently   provided   total   mean   squared   errors 
approximately   one   half  that    of   the   normal    distribution,    indi- 
cating a   superior    fit    in   the   right-hand    tail   area    of   the   dis- 
tribution.     Other   distributions   giving    good   fits    were   the 
Bernoulli-exponential    and   the  standard   exponential,    both 
yielding  total   mean   squared    errors   less    than   that    of    the   nor- 
mal   distribution. 


41 


B.  RECOMMENDATIONS 

The   results   of   this   study  indicate   the   Navy   should   con- 
tinue  to   use  the   Poisson  and   negative    binomial  distributions 
as   the   models   for  the    low-    and  medium-demand   classes   respec- 
tively.     Specifically,   the    Aviation   Supply   Office    should   re- 
evaluate  its   position   on   the   use   of   the   normal  distribution 
for    all   demand  classes  and    revert   to   the   Poisson    and   negative 
binomial  distributions  for    low-    and  medium-demand   classes   as 
before.      The  normal   distribution   tends    to   inflate    the    stock 
required   for   low-   and    medium-demand  items    for   a   specified 
level   of   protection   resulting  in   increased    safety    levels   and 
excessive  dollar   investment. 

The   Navy   should  consider   replacing   the   normal    distribu- 
tion   model    used    for    high-demand    items    with    the  compound 
Bernoulli-lcgnormal   distribution.      IE    the    Bernoulli-lognormal 
distribution   is   considered    computationally   too   difficult    for 
a  large    inventory   system,    the  Bernoulli-exponential   could   be 
used    in    its   place   with   little  loss  of   effectiveness. 

C.  ADDITIONAL    RESEARCH 

Follow  on  research  in  the  following  areas  may  improve  and 
expand  upon  the  results  presented  above: 

1.   Optimize  the  divisions  of  demand  categories  (i.e.  low, 
medium,  high)  with  the  possible  use  of  the  probability 
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of  a    non-zero   demand    as  an   element   in  the 
determination, 

2.  Implement   the    Bernoal li-lognormal    model   in    FMSO's   5A 
simulator   to    establish   the   effect    of   the  change   on   the 
entire   supply   system, 

3.  Investigate   the    implications  of    the   assumption    that    the 
demand   process    is   stationary, 

4.  Evaluate  other   compound  distributions,    such   as   the 
Bernoulli -log-logistic,    which   is    particularly   appealing 
analytically. 


U3 


APPENDIX    A 


GLOSS  ART   OF    TERMS 


Prograa   Related   Itea:      An   item  of  stock    whose   demand   can   be 
predicted   from   the   value   of    a   specific   Navy    program,    i.e., 
flying   hours,    steaming  hours,    etc. 

Reorder   Point:      The  on  hand    stock   quantity   that   when    reached, 
triggers   an   order   for   replenishment   of    stock   material.      The 
reorder   Point   is  the    expected   demand   during   procurement   lead- 
time   plus  the   safety    level. 

Risk:      Probability   of    a   stockout    during    leadtime. 

Safety   Level:      The   quantity    of   material    which   is    required    to 
be   on   hand    to   permit    continued  operation    in   the   event    of 
minor   interuptions   on    normal    replenishment   or   unpredictable 
fluctuations   in    demand.      The   safety   level   determined    is 
structured   so   as   to   minimize   time-weighted,    essentiality- 
weighted   requisitions    short. 

Stockout:      A   condition  that    exists   when   the   on   hand   inventory 
is   insufficient    to   fill   the    current   demand   requirements. 
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Systen   Perforaance:       A  subjective  measure   maximized   whar.   the 
time-weighted,    essentiality -weighted    requisitions    short    is 
minimized.      Time   weighting    is  the   consideration   of   the   aver- 
age  number    of    days    delay   in    the    availability   of   material, 
essentiality-weighting  is  the  consideration   of  the   relative 
essentiality  of   each    item.       Requisitions   short   are   requisi- 
tions   for   which    material    is    not    available. 
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APPENDIX_B 

gQNSgMAB|,E  ,ITEM.T  ASgiAT E D_R E SUL TS 
aO      REPETITIONS    FOR    EACH    OF         61    PSEDDO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBUTION:       POISSON 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^.80  O^^SS  0^90  0^95 

MEAN  0.750000      0.812500  0.854006  0.913924      0.956069 

STD    DEV  0.008766      0.007902  0.007148  0.005678       0.004149 

MSE  0.0  0.000156  0.000018  0.000290      0.000100 

TOTAL    MSE:  5.64E-04 

DISTRIBOTION:       NEGATIVE,.  BIN  OMIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0j^75  0^0  0^85  0^9  0  0^95 

MEAN  0.750000      0.812500  0.854006  0.913924       0.948994 

STD    DEV  0.008766      0.007902  0.007148  0.005678       0.004454 

MSE  0.0  0.000156  0.000018  0.000290      0.000322 

TOTAL    MSE:  7.86E-04 
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40       REPETITIONS    FOP.    EACH    OF         6  1    PSEUDO    SAMPLES 
LOR  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBUTION:       NORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0.75  0^80  0-.85  0^90  0^95 

MEAN  0.953221      0.953221  0.953221  0.953221      0.953221 

STD    DEV  0.004275      0.004275  0.004275  0.004275      0.004275 

MSB  0.042319      0.024488  0.011662  0.003835       0.001006 

TOTAL    MSE:  833.09E-04 

DISTRIBUTION:       EXPONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^80  0^85  0^90  0^95 

MEAN  0.953221       0.953221  0.953221  0.953221       0.953221 

STD    DEV  0.004275      0.004275  0.004275  0.004275       0.004275 

MSE  0.042319      0.024488  0.011662  0.003835       0.001006 

TOTAL    MSE:  833.09E-04 

DISTRIBUTION:       LOGNORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0j_80  0^85  0^19,  0^95 

MEAN  0.978404       0.995246  0.995246  1.000000       1.000000 

STD    DEV  0.002943       0.001392  0.001392  0.0  0.0 

MSE  0.052629      0.038292  0.021267  0.009998      0.002500 

TOTAL    MSE:       1246.85E-04 


47 


40       REPETITIONS    FOR    EACH    OF         6  1     PSSUDO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBUTION:       LAPLACE 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^25  2^8  0  0-.85  0^.90  0j^95 

MEAN  0.953221       0.953221  0.953221  0.953221       0.953221 

STD    DEV  0.004275      0.004275  0.004275  0.004275      0.004275 

MSE  0.042319      0.024488  0.011662  0.003835      0.001006 

TOTAL    MSE:  833.09E-04 

DISTRIBUTION:       LOGISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^.7  5  0^0  0^8  5  0^90  0.95 

MEAN  0.953221       0.953221  0.953221  0.953221       0.953221 

STD    DEV  0.004275      0.004275  0.004275  0.004275      0.004275 

MSE  0.042319      0.024488  0.011662  0.003835      0.001006 

TOTAL    MSE:  833.09E-04 

DISTRIBUTION:       BERNOULLI::EX  PONENTI AL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES;  0^75  0^80  0^85  g-.90  0 .  95 

MEAN  0.750000      0.812500       0.854006  0.913924       0.944939 

STD    DEV  0.008766      0.007902      0.007148  0.005673      0.004618 

MSE  0.0  0.000156       0.000018  0.000290      0.000613 

TOTAL    MSE:  10.77E-04 
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40       REPETITIONS    FOE    EACH    OF         61    PSEUDO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBUTION:       BERNOU  LLI~LO  GNORM  AL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^80  0^.85  0^.90  0^95 

MEAN  0.750000      0.312500       0.854006      0.913924       0.944939 

STD    DEV  0.008766      0.007902      0.007148      0.005678      0.004618 

MSE  0.0  0.000156       0.000018      0.000290      0.000613 

TOTAL    MSE:  10.77E-04 

DISTRIBUTION:       BERNOU  L  LI  ^OG I  STIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^5  0^80  0^.85  0-.90  0j^95 

MEAN  0.750000      0.812500      0.854006      0.913924  0.944939 

STD    DEV  0.008766      0.007902      0.007143      0.005678  0.004618 

MSE  0.0  0.000156       0.000018      0.000290  0.000613 

TOTAL    MSE:  10.77E-04 
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no       REPETITIONS    FOR    EACH    OF      753    PSEUDO    SAMPLES 
MEDIUM    DEMAND    RANGE:     1.0    <    D    <    20.0    PER    YEAR 

DISTRIBUTION:       PCISSON 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0-.75  0  .80  0-.85  0^90  0^95 

MEAN  0.748131       0.807719  0.843932  0.893923      0.925303 

STD    DE7  0.002501       0.002271  0.002091  0.001774      0.001515 

MSE  0.000138      0.000398  0.000603  0.001423      0.002266 

TOTAL    MSE:  48.28E-04 

DISTEIBOTION:       NEGATIVE    BINOMIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0-.75  2^80  2^.85  0.^10.  2i.21 

MEAN  0.746099      0.805430  0.843400  0.900210       0.943628 

STD    DEV  0.002508      0.002281  0.002094  0.001727      0.001329 

MSE  0.000332      0.000522  0.000589  0.000807      0.000606 

TOTAL    MSE:  28.55E-04 

DISTRIBDTION:       NORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0^90  0^95 

MEAN  0.888030      0.896699  0.907882  0.920152      0.935573 

STD    DEV  0.001817      0.001754  0.001666  0.001562      0.001415 

MSE  0.024771       0.014168  0.007327  0.003509       0.002400 

TOTAL    MSE:  521.74E-04 
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40       REPETITIONS    FOR    EACH    OF      753    PSEUDO    SAMPLES 
MEDIUM    DEMAND    RANGE:     1.0    <    D    <    20.0    PER    YEAR 

DISTRIBUTION:       EXPONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  O^JO  0^85  0^9  0  0^95 

MEAN  0.852981      0.861003  0.871591  0.884311       0.898352 

STD    DE7  0.002040      0.001993  0.001928  0.001843      0-001741 

MSB  0.020357      0-011779  0.006819  0.004953      0.006210 

TOTAL    MSE:  501.  19E-04 

DISTRIBUTION:       LOGNORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^.75  0^0  0^85  0^9  0  0-^5 

MEAN  0.953018      0.962552  0.970638  0.980623       0.990731 

STD    DEV  0.001219      0-001094  0.000973  0.000794      0.000552 

MSE  0.043846       0.028367  0.015977  0.007411       0.002040 

TOTAL    MSE:  976.41E-04 

DISTRIBUTION:       LAZLACE 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES;  0^75  0^0  0^85  0^9  0  0^95 

MEAN  0.884254       0.893767  0.904733  0.918110       0.936078 

STD    DEV  0.001843       0.001775  0.001692  0.001580      0.001409 

MSE  0.024079       0.013865  0.007191  0.003573       0,002363 

TOTAL    MSE:  510.  72E-04 
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t»0       REPETITIONS    FOR    EACH    OF      75  3    PSEUDO    SAMPLES 
MEDIUM    DEMAND    RANGE:      1.0    <    D    <    20.0    PER    YEAR 

DISTRIBUTION:       LOGISTIC 

THEORETICAL    PERCENTILES     (P) 

ESTIMATES:  Oj^Il  Q,j,M  ls.^1  0.^10,  2».21 

MEAN  0.875036      0.886687  0.896908  0.913037  0.934834 

STD    DEV  0.001905      0.001826  0.001752  0.001624  0.001422 

MSE  0.022871       0.013193  0.006955  0.003690  0.002447 

TOTAL   MSE:  491.56E-04 

DISTRIBUTION:       BERNODLLI^JX  PONENTIAI- 
THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^80  0^85  0^9  0  0^95 

MEAN  0.738072      0.795951       0.840430  0.896563      0.950685 

SID    DEV  0.002533      0.002322      0.002110  0.001755      0.001248 

MSE  0.002132      0.002338      0.002267  0.002224      0.001327 

TOTAL    MSE:  102.88E-04 

DISTRIBUTION:       3ERN0U LLI-LOGNORM AL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0^90  0^95 

MEAN  0.749344       0.809159       0.849391  0.903886       0.948263 

STD    DEV  0.002497      0.002264       0.002061  0.001698       0.001276 

MSE  0.001456      0.001751       0.001730  0.001866      0.001258 

TOTAL    MSE:  80.60E-04 
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40       REPETITIONS    FOR    EACH    OF      753    PSEODO    SAMPLES 
MEDIUM    DEMAND    RANGE:     1.0    <    D    <    20.0    PER    YEAR 

DISTRIBUTION:       BERNOU  LLI::LO  GISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  Oill  0^8  0  0^.8  5  0^9  0  0j^95 

MEAN  0.760006      0.821542       0.864176  0.916219      0.955416 

STD    DE7  0.002461      0.002206      0.001974  0.001596      0.001189 

MSE  0.002256      0.002499      0.002304  0.002030      0.001241 

TOTAL    MSE:  103.30E-04 
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40       REPETITIONS    FOR    EACH    OF      U29    PSEUDO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D      PER    YEAR 

DISTRIBUTION:       POISSON 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^7  5  Q^O  0.2.^1  ^s.2^  Q.^21 

MEAN  0.689018      0.730300  0.753320  0.780038      0.800697 

STD    DEV  0.003534      0.003388  0.003291  0.003162      0.003050 

MSE  0.024843      0.027146  0.033935  0.042154      0,051388 

TOTAL    MSE:       1794.66E-04 

DISTRIBUTION:       NEGATIVE    BINOMIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^8  0  0^85  0^90  0^95 

MEAN  0.724173       0.747130       0.791313  0.848041       0.896526 

STD    DEV  0.003412      0.003318      0.003102  0.002740       0.002325 

MSE  0.003047       0.024647      0.026316  0.029773      0.031637 

TOTAL    MSE:       1159,20E-04 

DISTRIBUTION:       NORMJIj 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0^.90  0^95 

MEAN  0.850264      0.871100  0.890154  0.910670       0.934248 

STD    DEV  0.002724      0.002558  0.002387  0.002177      0.001892 

MSE  0.016622      0.010312  0.005815  0.003155      0.002227 

TOTAL    MSE:         381.32E-04 
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UO       REPETITIONS    FOR    EACH    OF      U29    PSEODO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D       PEE    YEAR 

DISTRIBUTION:       EXPONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^8  0  2^85  0^9  0  0^,95 

MEAN  0.784610      0.814275  0.846685  0.881207      0-918702 

STD    DEV  0.003138      0.002969  0.002750  0.002470      0.002086 

MSB  0.009534      0.006885  0.005219  0.w04384       0.003800 

TOTAL    MSE:  298.21E-04 

DISTRIBUTION:       LOGNORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^80  0^8  5  0^9  0  0.95 

MEAN  0.855867      0.886448  0.918614  0.949144       0.973288 

STD    DEV  0.002681       0.002422  0.002087  0.001677      0.001113 

MSE  0.019028      0.013031  0.008113  0.004389      0.001639 

TOTAL    MSE:  462.00E-04 

DISTRIBUTION:       LAPLACE 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^5  0^80  0^85  0.90  0^95 

MEAN  0.843416      0.863600  0.885612  0,907473       0.934956 

STD    DEV  0.002774       0.002620  0.002430  0.002212       0.001883 

MSE  0.015775      0.009709  0.005680  0.003248       0.002167 

TOTAL    MSE:  365.79E-04 
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40       REPETITIONS    FOR    EACH    OF      429    PSEUDO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D      PER    YEAH 

DISTRIBUTION:       LOGISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0.75  0^.80  0^85  0^90  0^95 

MEAN  0.822132      0.846791  0.871709  0.900110      0.933755 

STD    DEV  0.002919      0.002750  0.002553  0.002289      0.001899 

MSE  0.013721       0.008841  0.005633  0.003523      0.002259 

TOTAL    MSE:  339.76E-04 

DISTRIBUTION:       BERNOULLI~EX  PONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^80  0^.85  0^.90  0^95 

MEAN  0.759754       0.816600       0.867492  0.914225       0.955101 

STD    DEV  0.003261       0.002954       0.002588  0.002138       0.001531 

MSE  0.005828      0.005034      0.003977  0.002618      0.001432 

TOTAL    MSE:  188.90E-04 

DISTRIBUTION:       BERNOO ILI-LQGNQRM  AL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^,75  0.^80  O^^SS  0^90  0^95 

MEAN  0.745078       0.797614       0.847634  0.901963      0.953670 

STD    DEV  0.003327      0.003067       0.002743  0.002270       0.001605 

MSE  0.004768      0.004403       0.003632  0.002613       0.001409 

TOTAL    MSE:  168.25E-04 
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40       REPETITIONS    FOR    EACH    OF      U29    FSEUDO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D      PER    YEAR 

DISTRIBUTION:       BERNOULLI -LOGISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIjIATES:  0^75  0^80  0^85  0^9  0  0;^95 

MEAN  0.828539      0.86U722      0.892571  0.921326      0  .9U7036 

STD    DEV  0.002877      0.002611       0.00236U  0.002055      0.001710 

MSE  0.012517      0.008862      0.005485  0.003020      0.001538 

TOTAL    MSE:         314.23E'04 
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APPJNDIX_C 

E2P AI RA BLE_I T EM-IiBOLil SD_R 2 SUITS 
40      REPETITIONS    FOR    EACH    OF      137    PSEODO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBOTION:       POISSON 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0^90  0.95 

MEAN  0.750000      0.812U96  0.353613  0.912594      0.955936 

STD    DEV  0.005849      0.005273  0.004775  0.003815      0.002772 

MSB  0.0  0.000156  0.000019  0.000311       0.000099 

TOTAL    MSE:  5.86E-04 

DISTRIBUTION:       NEGATIVE    BINOMIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0.^80  0^85  0^9  0  0^95 

MEAN  0.750000      0.812496  0.853613  0.912594      0.950753 

STD    DEV  0.005849      0.005273  0.004775  0.003815       0.002923 

MSE  0.0  0.000156  0.000019  0.000311       0.000270 

TOTAL    MSE:  7.57E-04 
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no       REPETITIONS    FOR    EACH    OF       137    PSEDDO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBUTION:       NORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0-.85  0-.22  2^11 

MEAN  0.949164      0.949164  0.949164  0.949164      0.949164 

STD    DEV  0.002967      0.002967  0.002967  0.002967      0.002967 

MSE  0.040918      0.023465  0.011016  0.003568      0.001118 

TOTAL    MSE:  800.84E-04 

DISTRIBUTION:    EXPONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^80  0^85  0^.90  0.95 

MEAN  0.949164      0.949164  0.949164  0.949-164      0.949164 

STD    DEV  0,002967      0.002967  0.002967  0.002967      0.002967 

MSE  0.040918      0.023465  0.011016  0.003568      0.001118 

TOTAL    MSE:  800.84E-04 

DISTRIBUTION:       LQGNORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^.75  0^.80  0^85  0^.90  0^95 

MEAN  0.978874      0.991028  0.991028  1.000000       1.000000 

STD    DEV  0.001943      0.001274  0.001274  0.0  0.0 

MSE  0.052955      0.036815  0.020201  0.009996      0.002500 

TOTAL    MSE:        1224.66E-04 
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40       REPETITIONS    FOR    EACH    OF       137    PSEODO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <   D    <       1.0    PER    YEAR 

DISTRIBUTION:       LAPLACE 

THEORETICAL    PERCENTILES     (P) 

ESTIMATES:                 0-.75                0^0                0^85  0^9  0  0^95 

MEAN                        0.949164      0.949164  0.949164  0-949164      0.949164 

STD    DEV                0.002967      0.002967  0.002967  0.002967      0.002967 

MSE                          0.040918      0.023465  0.011016  0.003568      0.001118 
TOTAL    MSE:          800.84E-04 

DISTRIBUTION:       LOGISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0.90  0^95 

MEAN  0.949164      0.949164  0.949164  0.949164       0.949164 

STD    DEV  0.002967      0.002967  0.002967  0.002967      0.002967 

MSE  0.040918      0.023465  0.011016  0.003563       0.001113 

TOTAL    MSE:         800.84E-04 

DISTRIBUTION:       BERNOULLI:iSX  PONENTI AL 

THEORETICAL    PERCENTILES     { P) 
ESTIMATES;  0^75  0^0  0^85  0;^90  0^95 

MEAN  0.750000      0.812496      0.853613  0.912594      0.942635 

STD    DEV  0.005849      0,005273      0.004775  0.003815      0.003141 

MSE  0.0  0.000156       0.000019  0.000311       0.000804 

TOTAL    MSE:  12. 91 E- 04 
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UO       REPETITIONS    FOR    EACH    OF       137    PSEDDO    SAMPLES 
LOW  DEMAND    RANGE:        0.0       <    D    <       1.0    PER    YEAR 

DISTRIBOTION:       BERNODLLI-LOGNORM AL 

THEORETICAL    PERCENTILES     (P) 

ESTIMATES:                0^75                0^0                2^85  0^9  0                0.95 

MEAN                       0.750000      0.812496      0.853613  0.91259a      0.942635 

STD    DEV                0.005849      0.005273      0.004775  0.003815      0.003141 

MSB                          0.0                  0.000156       0.000019  0.000311       0.000804 
TOTAL    MSE:             12.91E-04 

DISTRIBUTION:       BERNOU  LLI-LO  GIST IC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES;  0^75  0^8  0  0^85  0^90  0^.95 

MEAN  0.750000      0.812496       0.853613  0.912594      0.942635 

STD    DEV  0.005849      0.005273       0.004775  0.003815       0.003141 

MSE  0.0  0.000156      0.000019  0.000311       0.000804 

TOTAL    MSE:  12.91E-04 
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aO       REPETITIONS    FOR    EACH    OF    1487    PSEUDO    SAMPLES 
MEDIUM    DEMAND    RANGE:     1.0    <    D    <    20.0    PER    YEAR 

DISTRIBOTION:       POISSON 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  O^^O  0^,85  0^.90  0.95 

MEAN  0.747034      0.806323  0.843928  0.894732      0.929457 

STD    DEV  0.001782      0.001620  0.001488  0.001258      0.001050 

MSE  0.000238      0.000507  0.000582  0.001321       0.001836 

TOTAL    MSE:  44.84E-04 

DISTRIBOTION:       NEGATI  VE_BIN  OMIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^21  0^80  0^85  0^90  0^95 

MEAN  0.744688       0.304450  0.843869  0.900645       0.944953 

STD    DEV  0.001788      0.001626  0.001488  0.001227      0.000935 

MSE  0.000444      0.000598  0.000533  0.000788      0.000551 

TOTAL    MSE:  29.14E-04 

DISTRIBUTION:       NORMAL 

THEORETICAL    PERCENTILES     (P) 

ESTIMATES:                 0^7  5                0^0                0^85  0^9  0  0^95 

MEAN                        0.861474       0.875334  0.890316  0.905783       0.927155 

STD    DEV                0.001416      0.001354  0.001281  0.001198      0.001066 

MSE                           0.019444       0.011363  0.006084  0.003600       0.002944 
TOTAL    MSE:          434.35E-04 
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UO       REPETITIONS    FOR    EACH    OF    1487    PSSODO    SAMPLES 
MEDIUM    DEMAND    RANGE:     1,0    <    D    <    20.0    PER    YEAR 

DISTRIBUTION:    EXPONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES;  0^75  0^80  0-.85  9,^20  0^95 

MEAN  0.823217       0.834022  0.845365  0.864271       0.886389 

STD    DEV  0.001564      0.001526  0.001482  0.001404       0.001301 

MSE  0.016435       0.010140  0.007542  0.006691       0.007993 

TOTAL    MSE:         488.02E-04 

DISTRIBUTION:       LOGNORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATJS:  0^75  0-^0  1».3.1  0^90  0^5 

MEAN  0.946618      0.955351  0.964288  0-975138      0.987135 

STD    DEV  0.000922       0.000847  0.000761  0-000638       0.000462 

MSE  0.041761       0,026454  0.014803  0.006692       0.001890 

TOTAL    MSE:  916.00E-04 

DISTRIBUTION:       LAfiiCE 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0-8  0  0^85  0^90  0^95 

MEAN  0,856232      0.370161  0.885855  0.902481       0.928262 

STD    DEV  0.001439      0.001373  0.001304  0.001216       0.001058 

MSE  0,018823      0,011026  0,006113  0.003736      0.002851 

TOTAL    MSE:  425.48E-04 


63 


aO       REPETITIONS    FOR    EACH    OF    1487    PSEUDO    SAMPLES 
MEDIOM   DEMAND    RANGE:     1.0    <   D    <    20.0    PER    YEAR 

DISTRIBOTION:       LOGISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0.75  0j_80  0^85  0^9  0  0^95 

MEAN  0.845409      0.858922  0.875792  0.896303       0.925823 

STD    DE7  0.001482      0.001427  0.001352  0.001250      0.001075 

HSS  0.017690      0.010603  0.006269  0.004086       0.003093 

TOTAL    MSE:  417.42E-04 

DISTRIBUTION:       BERNOOLLI-EX  PONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0^90  0j^95 

MEAN  0.730065      0.737917       0.834508  0.893026      0.954267 

STD    DEV  0.001820      0.001676       0.001524  0.001267      0.000857 

MSE  0.003437      0.003608       0.003449  0.003021       0-001465 

TOTAL    MSE:  149.79E-04 

DISTRIBaTION:       BERNOUILIj^LOGNORM  AL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^.75  0^0  0^85  0.9  0  0^95 

MEAN  0.750179       0.808666       0.852134  0.904194       0.949477 

STD    DEV  0.001775      0.001613       0.001455  0.001207      0.000898 

MSE  0.002226      0.002681       0.002600  0.002327      0.001410 

TOTAL    MSE:  112.44E-04 
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40       REPETITIONS    FOR    EACH    OF    1487    PSEODO    SAMPLES 
MEDIUM    DEMAND    RANGE:     1.0    <    D    <    20.0    PER    YEAR 

DISTRIBUTION:       BERNOU  LLIiLO  gIST IC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^8  0  O2.85  0;^90  0^95 

MEAN  0.762406      0.822694       0.866844  0.915851       0.953988 

STD    DEV  0.001745      0.001566      0.001393  0.001138      0.000859 

MS2  0.003183      0.003523      0.003107  0.002468      0.001379 

TOTAL    MSE:  136.61E-04 
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aO       REPETITIONS    FOH    EACH    OF      356    PSEDDO    SAMPLES 
HIGH         DEMAND    RANGE;     20.0    <    D       PER    YEAS 

DISTRIBUTION:       POISSON 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^8  0  0i85  0^9  0  0^.95 

MEAN  0.695829      0.735554  0.762330  0.791447      0.821025 

STD    DEV  0.003855      0.003696  0-003567  0.003405      0.003212 

MSE  0.022812      0.026569  0-032141  0.038638      0.044655 

TOTAL    MSE:       1648.16E-04 

DISTRIBUTION:       NEGATIVE    BINOMIAL 

THEORETICAL    PERCENTILES     (P) 

ESTIMATES:                 0^75                0^.80                0^85  0^90  0^.95 

MEAN                       0.722861       0-755975  0.799852  0.857560      0.904930 

STD    DEV                0.003751       0.003599  0.003353  0.002929      0.002458 

MSE                          0.005285      0.024724  0.026795  0.028699      0.030113 
TOTAL    MSE:       1  156.  16E-04 

DISTRIBUTION:       NORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^7  5  0^.80  0^85  0^90  0^95 

MEAN  0-836334      0.859810  0.883941  0.905503       0.932035 

STD'  DEV  0.003100      0.002909  0.002684  0.002451       0.002109 

MSE  0.014669      0.009278  0.005429  0.003190      0.002319 

TOTAL    MSE:  348.85E-04 
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aO       REPETITIONS    FOR    EACH    OF      356    PSEUDO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D      PER    YEAR 

DISTRIBUTION:    EXPONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^.25  O^M  Qi^S  0^9  0  0^95 

MEAN  0.780940      0.818187  0.856974  0.898153      0.935255 

STD    DEV  0.003466      0.003232  0.002934  0.002535       0.002062 

MSB  0.009578      0.007296  0.005229  0.003830      0.002968 

TOTAL    MSE:  289.01E-04 

DISTRIBUTION:       LOGNORMAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  O^Jl  2^9.  2  5.8  5  0^9  0  0^95 

MEAN  0.819119      0.855337  0.895522  0.933281       0.971445 

STD    DEV  0.003226      0.002948  0.002563  0.002091       0.001396 

MSE  0.013003      0.008704  0.005985  0.003420      0.001454 

TOTAL    MSE:  325.65E-04 

DISTRIBUTION:       LAPLACE 

THEORETICAL    PERCENTILES     (?) 
ESTIMATES:  0^75  0^8  0  0^85  0^9  0  0^95 

MEAN  0.827652       0.852088  0.877334  0.902360       0.932870 

STD    DEV  0.003165       0.002975  0.002749  0.002487      0.002097 

MSE  0.013659      0.008798  0.005346  0.003302       0.002267 

TOTAL    MSE:  333.71E-04 
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no       REPETITIONS    FOR    EACH    OF      356    PSEUDO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D       PER    YEAR 

DISTEIBaTION:       LOGISTIC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^0  0^85  0^90  0^95 

MEAN  0.805812      0.833134  0.860882  0.892308      0.930748 

STD    DE7  0.003315      0.003125  0.002900  0.002598      0.002128 

MSE  0.012232      0.008347  0.005666  0.003803      0.002397 

TOTAL    MSE:  324-aUE-04 

DISTRIBUTION:       BERNOaLLI-EX PONENTIAL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^.25  0-JO  0^85  0.9  0  0^95 

MEAN  0.770942       0.829521       0.877192  0.921468       0.960362 

STD    DEV  0.003522      0.003151       0.002750  0.002254      0.001635 

MSE  0.007108      0,005916       0.004548  0.003092      0.001516 

TOTAL    MSE:  221.80E-04 

DISTRIBOTION:       BERNOU LLI-LOGNORM AL 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  0^75  0^30  0^8  5  0^9  0  0^95 

MEAN  0. 746108      0.797150       0.347197  0.901039      0.953617 

STD    DEV  0-003647      0.003370       0.003015  0.002502      0.001762 

MSE  0.005655       0.004764       0.003680  0.002809      0.001469 

TOTAL    MSE:  183.77E-04 
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UO       REPETITIONS    FOR    EACH    OF      356    PSEUDO    SAMPLES 
HIGH         DEMAND    RANGE:     20.0    <    D      PER    YEAR 

DISTRIBUTION:       BERNOUILI-LO  GIST IC 

THEORETICAL    PERCENTILES     (P) 
ESTIMATES:  O5.75  1^0.  l:.ll  Q.S.1Q.  9,£.ll 

MEAN  0.82983U      0.861163       0.886755  0.914062      0.941856 

STD    DEV  0.0031U9      0.002898      0.002656  0.002349      0.001961 

MSB  0.013074       0,008963      0.005495  0.003143       0.001722 

TOTAL    MSE:  323.98E-04 


69 


APPENDIX    D 

ESTIMATION  .CF    LOGNORMAL    MEA  N.,AND..VARI  ANCE_FROM  ,  UNTRANS  FORMED 

DATA 

The    usual    method   of   estimating   the    sample   mean   and   vari- 
ance   of   a   random   sample   assumed   to  come    from   a   lognormal   dis- 
tribution   is   to   take   the   log   transform   of   the   data   and   then 
compute   the  standard    method    of   moments   estimates.      This   is 
time   consuming   on   a   computer  as    the    logarithm   must    be   taken 
for    each  observation   in  the    sample.       Formulae   to    compute  the 
mean    and   variance   of   the   untransformed    data    from    the    mean   and 
variance   of  the   transformed    data   can    be    derived   for   the   log- 
normal  distribution. 

If   Y    is   distributed   nor  mal(Uy  ,ay)  ,    then   X    is    distributed 
lognormal   if   Y   =   ln(X).      The   mean   of    X    [Ref.    10]    is    expressed 
as: 

Uy   +  -~a^ 
^x   =    €  (20) 

and   the    variance   as: 

aj    =  e     ^  M  e^    -ij.  (21) 

Equations    (20)    and    (2  1)    are    the    reverse    of   that   required    in 
this   case   since    y^    and   a^   are  obtained    from   the   transformed 
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data.      What   is   needed   then,    is  to  solve    for  \iy   and    qs    in 
tens   of    y^  and    aj-      Substituting    (20)    into    (21)     results   in: 

which   when    solved   for    e      ,    taking   logs    and   simplifying,    gives 
the    solution    of    a^    as: 


,3  =  iJsLt/K 


(23) 


Substituting    (23)    back  into    the    log  of    (20)  ,    solving    for    Uy 
and   simplifying    gives   the   solution   of    ijy   as: 


(24) 


Thus   the   mean   and   variance    cf   the   transformed   data   may   be 
estimated   from   the   sample   mean   and   variance   of  the   untrans- 
formed   data   reducing   significantly   the   computational 
requirements. 
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